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20.1 SUMMARY 


Ecotoxicity of chemicals at the population level may be assessed on the basis of 
experiments in which test organisms are exposed for the whole length of their life-cycle, 
and observations are made for growth, development, reproduction and survival. The 
data obtained from these experiments may be analysed using demographic techniques 
that allow population parameters, such as the intrinsic rate of increase and the biomass 
turnover rate, to be estimated. Ag--structure and stage-structure models, using either 
continuous time or discrete time, may be used for the analysis of ecotoxicological data. 
This chapter evaluates the application of demographic approaches in the ecotoxicologi- 
cal literature; a total of 29 studies are reviewed. Short-lived aquatic invertebrates are the 
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most popular test organisms. However, more recently, soil invertebrates have been 
increasingly used. The studies show that, in contrast to some suggestions in the literature, 
population parameters are not necessarily more sensitive than individual parameters and 
the most sensitive individual endpoint may not be the most relevant in terms of 
population level effects. The use of population parameters lies in their ecological rel- 
evance, and the possibility of summarizing a variety of effects in the life-cycle using 
a single endpoint. An "elasticity index of toxicant exposure" is proposed to measure the 
change of a life-cycle variable, and its effects on fitness, with a change of the exposure 
concentration. Although theory shows that the tr: de-off structure among life-cycle 
variables is important, empirical studies have not been able to reveal these trade-offs 
clearly. The application of demographic methods in ecotoxicology may contribute to 
improving the ecological relevance of ecotoxicological risk assessment. 


20.2 INTRODUCTION 


A complete impression of the ecotoxicity of a chemical to a single species is 
obtained from life-cycle experiments. In these experiments individuals or groups 
of individuals are exposed to a series of constant concentrations, and observa- 
tions are made on various life-history events, such as (depending on the species) 
hatching, germination, growth, moulting, flowering, seed formation, egg laying, 
development, survival. Each of these observations is a potential endpoint in 
which the toxicity of the chemical may be expressed. Often one will find that the 
endpoints differ considerably in their response to toxicants. This leaves the risk 
assessor with the question: which criterion should be chosen for the protection 
of the species? 

The dilemma of sensitivity is illustrated in a theoretical way in Figure 20.1a. 
Three concentration-response relationships are shown, for different life-history 
traits, each with its own EC; (50% effect concentration) and NEL (no effect 
level). Usually survival is a less sensitive endpoint than body-growth or repro- 
duction, However, this depends on the mode of action of the chemical and the 
physiological response of the species. To protect the species, should one take the 
most sensitive criterion (in this case NEL for body-growth), or is it better to 
focus on one of the criteria, e.g. reproduction? A similar problem may arise when 
two species are compared (Figure 20.1b). If the two species differ in their basic 
reproductive capacity, but have the same EC;,, will the toxicant affect the two 
species in a similar manner? 

Population ecology, specifically life-history theory, offers a straightforward 
answer to the questions posed above: one should take the criterion that has the 
greatest influence on the intrinsic rate of population increase, or even better: one 
should integrate all life-history traits into this single parameter. Various 
ecotoxicological studies have followed this approach, often referred to as “life- 
table evaluation". This chapter will review these studies, discuss some related 
approaches, and evaluate their contribution to the improvement of ecotoxico- 
logical risk assessment. 
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Figure 20.1 Theoretical concentration-response relationships for one species with three endpoints 
differing in sensitivity (a), and two species (A and B) with the same toxicological sensitivity but 
a difference in control reproduction (b). NEL = no-effect level, LCsy = median lethal concentration, 
ECs» = 50% effect concentration. 


The question may be raised as to whether the study of population parameters 
might provide endpoints that are more sensitive than individual parameters 
(Suter and Donker 1993). This argument was initially formulated by Halbach et 
al. (1983), who called it the “looking glass effect": small, hardly perceptible, 
effects on individuals are magnified by the great number of individuals constitut- 
ing a population, and so become more easily visible. Recent literature, however, 
seems to agree on the fact that population parameters cannot be more sensitive 
than individual parameters, and, in fact, are sometimes less sensitive (Day and 
Kaushik 1987; Meyer et al. 1987; Barbour et al. 1989; Crommentuijn et al. 1993). 
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Figure 20.1 (Continued) 


Hence, the greatest value of the use of population parameters does not lie in 
their assumed sensitivity, but in their ecological relevance, and the possibility of 


summarizing the variety of possible effects in the course of a life-cycle by a single 
measure. 


20.3 DEMOGRAPHIC TECHNIQUES 


The classic demographic methodology starts with a summary of the life-history 
using two functions, I(x) and m(x), that define survival (fraction still alive, 
relative to birth), and fertility (no. offspring per time unit per parent) as 
a function of age, x. From these two functions, the Malthusian parameter, r, may 
be calculated by means of Euler's equation: 


i l(x)m(x)e7* dx = 1 (20.1) 


x-0 
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This equation relates to the basic demographic theorem (see, e.g., Keyfitz 
1968) that may be formulated as follows: 


Any age-structured population with time-invariant l(x) and m(x) schedules 
will approach a stable age distribution, independent of the initial conditions; 
the population will then increase at an exponential rate r (per time unit). 


From an ecological point of view, the relevance of equation 20.1 is contained 
in Fishers fundamental theorem of natural selection, which states that, in 
a population composed of genotypes with different life-histories, the types with 
the highest Malthusian parameter (r) tend to displace others from the popula- 
tion, and that the rate of change of the population growth rate is approximately 
equal to the additive genetic variance of r (Charlesworth 1980). In other words, 
ris a measure of fitness (Kozlowski 1993). This theorem forms the basis for the 
study of life-history evolution, in which an analysis is made of the changes in the 
life-history that maximize fitness, given certain trade-off patterns between life- 
history parameters (Stearns 1992). In an ecotoxicological context, the relevance 
of equation 20.1 may be summarized in the following statement: 


The population consequences of the effects of a toxicant on survival and 
fertility must be judged by the extent to which a reduction in survival or 
fertility affects the Malthusian parameter of the life-cycle. 


Despite its obvious ecological relevance, the use of r has a drawback in that it 
only measures the numerical aspects of the population (size in numbers, mortal- 
ity and fertility in numbers per time), not its qualitative aspects. At the individual 
level body-weight is another common endpoint, and its analogue at the popula- 
tion level is total biomass. Van Straalen (1985a) and Aldenberg (1986) discussed 
the concepts of biomass and productivity, as defined for populations structured 
by body-size instead of age. Following from this theory, Van Straalen and De 
Goede (1987) proposed the intrinsic biomass turnover rate, P/B, as a new 
population performance index that includes possible effects of a toxicant on 
body-growth. The intrinsic biomass turnover measures the productivity of the 
population, that is the rate at which biomass is produced, relative to the biomass 
present. Production of biomass may be a more relevant measure for the 
performance of a natural population than production of numbers, especially in 
relation to its functioning in trophic chains. The parameter P/B may be cal- 
culated from the age-specific survival rate, I(x), the age-specific fertility, m(x), and 
the age-specific body weight, w(x), by means of the following equation: 


i — dl je dx 
de Bx y (202) 


i l(x)w(x)e^" dx 
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The derivation of this equation, given in Van Straalen and De Goede (1987), is 
based on what may be called the weight-extended demographic theorem: 


Any structured population with time-invariant survival, fertility and body- 
growth will approach a stable weight distribution, independent of the initial 
conditions. The total biomass of the population (B), and its total rate of 
production (P) will then both increase exponentially at a rate r, while their 
quotient (P/B) will approach a constant. 


The calculation of r from real data arranged in a life-table is explained in 
Stearns (1992). It requires the use of a (simple) computer programme, because of 
the implicit expression of r in equation 20.1. In a similar manner, P/B may be 
calculated from a life-table supplemented with body-growth data (Van Straalen 
1985a; Van Straalen and De Goede 1987). 

The reference to "time invariance" in the above statements may be considered 
a serious limitation of the application of the theory, because the vital statistics of 
a population often show considerable variation, in response to changing envi- 
ronmental conditions. The value of the statements should therefore not be 
judged by the extent to which r and P/B are able to provide estimates for the real 
rates of increase and productivity in a population under field conditions. Rather, 
they must be considered instantaneous values that summarize the fitness of 
a species in a particular environment and indicate the potential for population 
growth in the future, if the present conditions persist. So it is more appropriate 
to interpret r and P/B as instantaneous performance indices, rather than as real 
predictions of future behavior. 

When temperature is the only variable affecting the changes in vital rates, the 
time-invariant approach may still be used after a suitable transformation of the 
time-scale. Van Straalen (1983) argued that, under certain conditions, any 
process whose rate changes with time due to temperature changes, can be made 
to proceed at a constant rate by warping the time-scale. The new time-scale, 
referred to as "physiological time", can be defined in terms of the temperature 
responses of the specific organism and the change in temperature with time. This 
approach has proven useful when modelling population development of small 
ectotherms, for which temperature is the dominating factor determining life- 
history events (Van Straalen 1985b; Stamou et al. 1993; Diekkrüger and Róske 
1995; Axelsen 1997). 

The mathematical theory underlying the derivation of equations 20.1 and 20.2 
considers time and age as continuous variables. The techniques used are those of 
calculus algebra. Alternatively, age may be considered a discrete variable. All 
events are then assumed to take place at certain moments in time. At regular 
intervals (of length T ) individuals are assumed to die, produce offspring, or 
jump from one class to the next, while nothing happens in between. The 
mathematical techniques used in this approach are those of matrix algebra. 
Development of the population may be described by a matrix equation in which 
the age distribution at a certain moment (considered a column vector) is 
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multiplied by a matrix L, to obtain the age distribution one time step later. The 
population projection matrix L is usually referred to as Leslie's matrix, and has 
a typical configuration, as follows: 


0 f f h 
p, 0 0 0 
EEUU 0 
t=l bm x (20.3) 
0 0 0 Pr-1 0 | 


All the elements of the projection matrix are zero, except for the off-diagonal 
elements p;, which represent the probabilities of survival from age-class i to age 
class i+ 1, and some of the elements in the first row, fi, which represent the 
number of offspring produced at age i and surviving to the next time step, when 
they enter the zeroth age class. Demographic theory shows that the projection 
matrix has one real dominant eigenvalue, usually denoted by 4. This parameter 
represents the (dimensionless) factor by which the population is multiplied per 
time interval T, when it is in a stable age distribution and growing exponen- 
tially. The parameter 4 is also referred to as the "finite rate of increase" and it is 
related to r according to r T — In 4, where T is the time step considered (usually 
unit of time). 

The matrix description of population development has the advantage that it is 
more flexible towards incorporating time-dependent processes. Projection 
matrix elements may be given different values for different time steps (e.g., a stop 
on reproduction in the winter period). An annual projection matrix may be 
composed of twelve monthly matrices, each describing the life-history traits in 
that month. This approach is therefore more suitable for predicting population 
changes under field conditions. The matrix approach is easy to formulate in 
computer programmes and is often used in simulation models. 

The Leslie matrix approach may be extended to include stage-structured 
populations rather than populations structured by age only. Stages may be real 
biological stages, such as the instars of arthropods, or they may be arbitrary 
classes of body-weight. Lefkovitch (1965), Vandermeer (1975), Longstaff (1977) 
and Van Groenendael et al. (1988) considered projection matrices for stage- 
structured populations. This theory is briefly summarized here. 

Compared to the classical Leslie approach, the matrix description of stage- 
structure changes poses a problem because the stages are of unequal duration 
and the individuals in the population do not all jump from one stage to another 
after each time interval. After each time step T there are three possibilities for an 
individual in stage i: to die, to stay in stage i or to jump to stage i + 1. Transition 
rates are denoted by qj; (probability to stay in stage i) and q;;. , (probability 
to jump to stage i -- 1). The projection matrix, M, then has the following 
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configuration: 


[ doo fh h h 7] 
doi: di). O 0 
O qui 42.2 0 
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M 0 © da 0 (20) 
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In contrast to the Leslie matrix, the survival probabilities per time interval 
T are now not equal to the off-diagonal elements q;,;+; but to qii + diis 
Schobben and Van Straalen (1987) showed that the matrix elements of the 


stage-structure projection matrix may be given the following biological inter- 
pretation: 


= pi(di -— y 
di 


(20.5) 


[E] 


where p, is the survival rate per time interval T, as above, and d; is the duration 
of stage i. The Leslie matrix L may be considered as a special case of the 
stage-structure projection matrix M. It can be seen from equations 20.5 that if all 
stage durations d; are the same and equal to the time interval T, the probability 
of staying in a stage, q,;, equals zero, qi i+; equals p; and the two matrices are 
identical. So, the description of stage-structured populations does not necessar- 
ily require the completely filled projection matrices proposed by Lefkovitch 
(1965). If it is assumed that individuals do not skip stages (that is, do not jump 
from stage i to stage i + 2 during time interval T), only one extra diagonal of 
elements is necessary to extend the Leslie matrix to stage-structured popula- 
tions. The prohibition to skip stages can easily be met by making the time 
interval T shorter than the shortest stage duration. 

In addition to the parameters r, P/B, and A, other population character- 
istics have been used in ecotoxicological studies for describing certain aspects 
of population change. Gentile et al. (1982) proposed the use of reproductive 
value, v(X), as a measure of the expected contribution of a female to future 
population growth. Various authors have used the carrying capacity, K, of the 
population, usually estimated after fitting the logistic growth model to experi- 
mental data. Halbach et al. (1983) introduced two other parameters, the fre- 
quency of oscillation (f), and the “pregnancy”, p, of population fluctuations. 
The latter index is a measure of the persistence of oscillations when these are 
analysed using autocorrelation functions. The indices K, f, and p do not have 


à strong theoretical underpinning, they are mere descriptors of population 
fluctuations. 
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20.4 POPULATION PERFORMANCE INDICES 
IN ECOTOXICOLOGICAL STUDIES 


The demographic techniques reviewed above and elsewhere (Caswell 1996; 
Kareiva et al. 1996) demonstrate that there is a good theoretical framework for 
the application of the life-history approach in ecotoxicology. To exemplify the 
approach, some examples from the literature will be reviewed in this section. 

Demographic analysis of ecotoxicity at the population level has been applied 
to a variety of different species. Table 20.1 lists studies reported in the literature, 
with a specification of the species and chemicals investigated. Among the 
species, freshwater invertebrates with short life-cycles are represented well, 
presumably because these species can easily be cultured and observations on 
life-history traits can be made within a relatively short experimental period. The 
most widely used population performance index is r, the intrinsic rate of 
population increase, derived from Equation 20.1, or estimated directly from 
populations growing exponentially in laboratory cultures. More recently, terres- 
trial ecotoxicology studies have also considered the demographic approach. 

By way of example, Figure 20.2 demonstrates the use of the intrinsic rate of 
population increase for populations of rotifers, cladocerans and algae, exposed 
to various concentrations of 3,4-dichloroaniline (Seitz and Ratte 1991; Ratte 
et al. 1992). The response curves clearly show the differential effects of the 
toxicant on the species investigated. The use of the intrinsic rate of population 
increase in this case made it possible to express toxicity in terms of an endpoint 
that is relevant to all five species. Although the toxicant may affect different 
life-history traits in different species, the effects can be compared on a common 
basis by using r. . 

Another lesson to be learnt from the studies listed in Table 20.1 is that the 
population effects of a toxicant, particularly when they reduce juvenile survival, 
greatly depend on food availability. This was shown in a study by Kooijman 
(1985). The toxicity of vanadate (affecting survival) to the rotifer Brachionus 
rubens, expressed in terms of r, was hardly visible at high food concentrations, 
but became clearly apparent at low food availability. Toxicity of 3,4-dichloro- 
aniline (affecting reproduction) was, however, independent of food supply. These 
observations are in line with predictions from a dynamic energy budget model 
(Kooijman and Metz 1984; Kooijman 1993). . 

An example of the matrix approach in ecotoxicology can be found in Schob- 
ben and Van Straalen (1987). In this study, a model was constructed for the 
population dynamics of the oribatid mite Platynothrus peltifer, which develops 
in six discrete stages: egg, larva, three nymphal stages, and the adult stage. Based 
on life-history information obtained from field studies, survival probabilities, 
stage durations, and reproduction activities were specified for each stage, and 
were assumed to vary per month. The annual population projection matrix was 
obtained by subsequently multiplying the monthly projection matrices. The 
eigenvalue of the annual projection matrix (4) represents the factor by which the 


Table 20.1 A historically ordered listing of studies in which demographic techniques were applied to analyse the population effects of toxicants, 


specifying species, chemicals, and population endpoints considered 


Reference Species Group System Chemical Population 
endpoint" 
Winner et al. (1977) Daphnia magna Cladocera freshwater copper r 
Marshall (1978a, b) Daphnia galeata Cladocera freshwater cadmium r, K 
Daniels and Allan (1981) Eurytemora affinis Copepoda marine dieldrin r 
Daphnia pulex Cladocera freshwater dieldrin r 
Allan and Daniels (1982) Eurytemora affinis Copepoda marine kepone r 
Gentile et al. (1982) Mysidopsis bahia Mysidacea marine mercury, nickel r, v(x) 
Halbach et al. (1983, 1984) Brachionus rubens Rotifera freshwater pentachlorophenol r, K, fp 
Kooijman (1985) Brachionus rubens Rotifera freshwater 3,4-dichloroaniline, r 
vanadate 
Bengtsson et al. (1985) Onychiurus armatus Collembola soil copper, lead r 
Van Leeuwen et al. (1985) Daphnia magna Cladocera freshwater cadmium r, K 
Van Leeuwen et al. (1986) Daphnia magna Cladocera freshwater bromide r, K 
Van Leeuwen et al. (1987) Daphnia magna Cladocera freshwater cadmium, mane., rK 
dichromate, vanadate, 
2,4-dichloroaniline, 
pentachlorophenol, 
TPBS, bromide, 
pentachlorobenzene 
Van Straalen and De Goede (1987) — Orchesella cincta Collembola soil cadmium P/B 
Day and Kaushik (1987) Daphnia galeata Cladocera freshwater fenvalerate È 
Meyer et al. (1987) Daphnia pulex Cladocera freshwater cadmium, copper r 
Schobben and Van Straalen (1987) — Platynothrus peltifer Oribatida soil cadmium À 
Barbour et al. (1989) Ceriodaphnia Cladocera freshwater complex effluents r 
Van Straalen et al. (1989) Orchesella cincta, Collembola soil cadmium r 
Platynothrus peltifer Oribatida 
Table 20.1 (Continued) 
Reference Species Group System Chemical Population 
endpoint ™ 
Van der Hoeven (1990) Daphnia magna Cladocera freshwater 3,4-dichloroaniline, r 
metavanadate 
Wong and Wong (1990) Moina macrocopa Cladocera freshwater cadmium r 
Seitz and Ratte (1991) Daphnia, Chydorus, Cladocera freshwater 3,4-dichloroaniline r 
Ceriodaphnia, Cladocera 
Keratella Rotifera 
Crommentuijn et al. (1993) Folsomia candida Collembola soil cadmium r 
Enserink et al. (1993) Daphnia magna Cladocera freshwater cadmium, chromium r 
Bechman (1994) Tisbe furcata Copepoda marine copper r 
Klüttgen and Ratte (1994) Daphnia magna Cladocera freshwater cadmium r 
Kammenga et al. (1996) Plectus acuminatus Nematoda soil cadmium r 
Klok and de Roos (1996), Lumbricus rubellus Lumbricidae soil copper r 
Klok et al. (1997) 
Hendriks and Enserink (1996) Phalacrocorax carbo Aves freshwater PCB-153 r 
Crommentujn et al. (1997) Folsomia candida Collembola soil cadmium, chlorpyrifos, r 
triphenyltin 
Laskowski (1997) Helix aspersa, Gastropoda terrestrial cadmium, copper, zinc, lead — 4 
Lithobius mutabilis Chilopoda 


Vr = intrinsic rate of increase, 4 = population multiplication factor ( = e'), P/B = intrinsic rate of biomass turnover, v(x) = reproductive value, K = carrying 


capacity, f = oscillation frequency, p = "pregnancy" of population fluctuations with time. See text for explanation. 
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Figure 20.2 Concentration-response relationships for the effect of 3,4-dichloroaniline (DCA) on the 
population growth rate (r) of five species of aquatic organisms. For each species, r is expressed as 
a percentage of the value in the blank treatment. After Seitz and Ratte (1991). 
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Figure 20.3 Calculated population increase, expressed as the annual replacement factor, A, of the 
oribatid mite Plat ynothrus peltifer, exposed to cadmium-polluted food. The factor À was estimated as 
the eigenvalue of the annual projection matrix (obtained after subsequently multiplying 12 monthly 


projection matrices), in which the life-cycle in the field was combined with laboratory observations 
on cadmium toxicity. After Schobben and Van Straalen (1987). 


population is multiplied yearly. To maintain a stable population over the years, 
À should equal unity. 


Figure 20.3 illustrates some of the results obtained when the matrix model 
for P. peltifer was combined with observations for the effects of cadmium on 
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reproduction, stage duration and survival. It may be seen that — E 
a rather drastic effect on this species, A is reduced below unity already by i e 
third exposure concentration. At this and all higher concentrations the acd 
tion will become extinct (A < 1 implies r « 0). The use of the matrix model in : is 
case provided a way to summarize a variety of complicated bae. al : 
nomena (development through stages of unequal duration, him a e 
ity, etc.) as a single measure for the performance of the population. t also 
illustrated that the effect of cadmium on 4 was more easily visible, compared to 
the effects on individual reproduction (cf. Van Straalen et al. 1989). 


20.5 SENSITIVITY ANALYSIS BASED ON LIFE-HISTORY 
RESPONSES TO TOXICANTS 


The concept of sensitivity is not an easy one. Meyer et al. (1987) net 
various aspects relating to the question of how to decide whether one variable 
responds more sensitively to a toxicant than another variable. From a toxico- 
logical point of view, sensitivity may be discussed in terms of the — 
sure-response graph. An endpoint is considered to be more sensitive il 
toxicological effect criterion (e.g, ECso) is positioned more to the left on t e 
exposure axis. In practice, however, sensitivity is also determined by the possi- 
bility of demonstrating an effect emerging above the noise in the data. It bh 
be possible, in principle, to have a variable that is very sensitive in a toxicolo- 
gical sense, but that has such a high degree of statistical uncertainty that it 
requires relatively high exposure levels to demonstrate a response ag 
significantly from the control. This is particularly relevant in the case o 
secondary parameters, such as the intrinsic rate of population increase, which 
are usually estimated from different sources of information, each contributing to 
increase the error. The statistical aspects of estimating confidence intervals for r, 
using numerical techniques, were discussed by Meyer et al. (1986). Although we 
recognize the practical importance of the statistical concept of sensitivity, it will 
be ignored for the purpose of this paper and sensitivity will be discussed in the 
classical toxicological sense. N 

When the life-cycle of an organism is specified in demographic terms, and 
equation 20.1, or its discrete time analog, is used to estimate ror 4, it becomes 
possible to analyse how the value of r (or 4) changes with changes in life-history 
traits. In an ecological context, this analysis simulates the effects of gene 
substitutions that affect a certain life-history trait. The response of r tells us how 
strong the selection pressure is on that particular trait. This type of sensitivity 
analysis dates back to Cole (1954), and has been elaborated intensively in the 
context of life-history theory (Charlesworth 1980; Stearns 1992; Caswell 1996). 
Caswell (1978) approached the problem by considering the relative change in 
À brought about by a small relative change in one of the matrix elements 
of the population projection matrix, if all the other matrix elements were held 
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constant. In mathematical terms, this means that one can consider the partial 
derivative of In A with respect to the natural logarithm of the matrix element as 
a measure of the "elasticity" of that matrix element (see also Van Groenendael et 
al. 1988). Elasticity of the a,; element is denoted as e;; and is defined as follows: 


dln 4 


= ain ai; (20.6) 


ej 
This type of analysis has led to the conclusion that fitness is sensitive to different 
life-history traits under different rates of population growth. In a rapidly 
expanding population, there is a strong selection pressure on age at maturity, 
while in a stable population, survival and fertility rates are more important 
variables (Stearns 1992). 

The sensitivity analysis from life-history theory was applied in an ecotoxico- 
logical context by Kammenga (1995) and Kammenga et al. (1996). Figure 20.4 
shows a theoretical illustration of the argument. In the right part of the figure, 
a normal concentration-response relationship is seen, for two different species, 
1 and 2. The effects of the toxicant are expressed in terms of some life-history 
trait (e.g., reproduction, survival, development). At a certain concentration of the 
toxicant, species 2 is affected to a greater extent than species 1. In the left part of 
the curve, graphs are shown that relate the intrinsic rate of population increase, 
r, to the life-history trait considered. The theoretical example is chosen such that 
ther value of species 1 reacts much more strongly to a decrease in the life-history 
trait than the r value of species 2. The result is that r of species 1 is affected to a 
greater extent than r of species 2, even though the two species are exposed to the 
same exposure concentration and species 2 is more sensitive with respect to the 
life-history trait considered. 
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Figure 20.4 Theoretical response curves for the effect of a toxicant on a life-history parameter (right) 
and the influence of this parameter on the intrinsic rate of population increase (left). Curves are 
shown for two species, | and 2, that are assumed to differ in sensitivity. 
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A sensitivity analysis of the type outlined above was applied by Kammenga 
et al. (1996) in experiments using the soil nematode Plectus acuminatus exposed 
to cadmium. Their analysis challenged the generally held opinion that popula- 
tion level effects of toxicants are determined by the effects on the most sensitive 
life-history trait. This is, for example, the rationale behind the development of 
the early life-stage test for fish (McKim 1977). However, Kammenga et al. (1995) 
showed that in the nematode experiments the effect of cadmium on r was not 
determined by the effect on the most sensitive life-history component, the 
reproductive period. In fact, a reduction of the reproductive period by 45% did 
not reduce fitness to any extent, while a prolongation of the juvenile period, the 
least sensitive trait for cadmium toxicity, had a clear effect on r. 

Similar analyses of sensitivity to toxicants have been applied to bird popula- 
tions, using demographic models (Young 1968; Emlen and Pikitch 1989; Meyer 
and Boyce 1994). These studies have emphasized the importance of different 
life-history traits to population growth rates, depending on the life-time of the 
bird species. In addition, Meyer and Boyce (1994) showed that to evaluate the 
effect of pesticides on avian populations, it is essential to estimate the variability 
in survival and fertility. 

The fitness effects of toxicants may be analysed by means of elasticity indices 
that measure not only the change of fitness with a small change of a life-history 
variable, as in equation 20.6, but also the response of fitness to a change in the 
exposure concentration (Kammenga 1995). In analogy to equation 20.6, one 
may define a new elasticity index, “toxicant exposure elasticity", as the partial 
derivative of In 2 with respect to the exposure concentration. Because concentra- 
tions act on a multiplicative scale rather than on a linear scale it is logical to 
consider a change in the logarithm of the concentration (In c), and define the 
toxicant exposure elasticity index, k, as: 


,  Olnà Ana, 
~ nai; Olnc 


(20.7) 


This index would be a direct measure of the change in fitness that is experienced 
by a species when the concentration of a toxicant in the environment increases 
by a small amount, relative to the concentration present (ln c = (1/c)éc). 
Another approach to sensitivity analysis was developed by Crommentuijn 
et al. (1995). Following Daniels and Allan (1981), Crommentuijn et al. (1995) 
considered the concentration-response relationships of survival and sublethal 
characteristics, such as reproduction, in a single graph (Figure 20.5). The ratio 
between the LCs, and the threshold concentration for reproduction was defined 
as the "sublethal sensitivity index", abbreviated SSI. A large value for this index 
indicates that reproduction is inhibited at a level far below the LC; (a type III 
response in the terminology of Daniels and Allan 1981), while a low value 
indicates that the organism maintains reproduction until it dies (type | response 
according to Daniels and Allan 1981). An inventory of SSI-values for several 
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Figure 20.5 Showing hypothetical concentration-response relationships for survival and reproduc- 
tion of the same species, from which the sublethal sensitivity index, SSI, according to Crommentuijn 
et al. (1995), may be derived as the quotient of LCso and the no-effect concentration (NEC) for 
reproduction. 


invertebrates exposed to cadmium showed that it varied between 280 and 3.1 
(Crommentuijn et al. 1995). SSI may be considered a measure of an organism's 
priorities when, under toxicant stress, it is forced to allocate energy either to 
reproduction or survival. It can be an interesting tool for the analysis of 
species-specific responses to toxicants, as related to the life-history strategy. 
Previous analyses have taken the r and K selection theory as a basis for 
classifying life-history related responses to toxicants, but with limited success 
(Neuhold 1987). 


20.6 TRADE-OFF RELATIONSHIPS AMONG 
DEMOGRAPHIC CHARACTERS 


In the previous section, the question of sensitivity was approached on a charac- 
ter-by-character basis, that is, the influence of one life-history trait on r or A was 
considered under the condition that all other traits remained the same. This is 
a rather unrealistic situation, because it is well known that a change in one 
variable will often evoke a correlated response in another variable. When the 
secondary response can be considered as (partly) compensatory to the first, the 
relationship between the two variables is called a trade-off-relationship. For 
example, a reduction in fertility may increase the probability of adult survival. In 
these cases, the question of sensitivity cannot be solved for one isolated charac- 
ter, but the life-history must be seen as a whole. 
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Figure 20.6 Correlation between cocoon production and body growth in experiments using the 
earthworm Eisenia andrei, cultured for 3 weeks in OECD artificial soil substrate. The data are from 
18 separate experiments (all without toxicants). After Van Gestel et al, (1992). 


As an example of a trade-off mechanism, Figure 20.6 gives the results of 
a series of experiments using the earthworm Eisenia andrei in an artificial soil 
substrate (Van Gestel et al. 1992). The production of cocoons in the controls of 
ecotoxicity experiments was plotted against the increase in body mass in the 
same group of animals. Figure 20.6 clearly shows that there is a negative 
correlation between these two variables: in experiments where the worms 
produced many cocoons, their body growth was slow, sometimes even negative. 
Conversely, toxicants that inhibit reproduction may increase the weight gain of 
these animals. A similar trade-off mechanism was noted by Crommentuijn et al. 
(1997) in experiments where the collembolan Folsomia candida was exposed to 
cadmium and triphenyltin. Both toxicants reduced body growth significantly, 
with a correlated reduction of egg output, but adult survival time increased with 
a decrease in growth. 

The trade-off structure among life-cycle variables was the basis for a theory 
formulated by Sibly and Calow (1989). This theory attempts to link the physi- 
ology of an organism with its life-history, and consequently with its population 
growth. Calow and Sibly (1990) visualized the theory by defining the state ofan 
organism in a “physiological state space” diagram. Production and metabolism 
were considered the most important variables for characterizing the physiol- 
ogical state (Figure 20.7). In the diagram, isoclines were drawn to indicate all 
combinations of production and metabolism that have equal values for the 
intrinsic rate of increase, r. Three different lay-outs of the r isoclines were 
proposed, corresponding to a stepwise dependence of production on metabol- 
ism (left), a linear decrease (middle) and a curvilinear decrease (right). When a 
toxicant moves the physiological state of an organism in a certain direction, the 
consequences for the population are determined by the trajectory it takes in the 
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Figure 20.7 Physiological state space diagrams, showing the theoretical relationship between 
metabolism and production, according to Calow and Sibly (1990). The isoclines indicate all 
combinations of metabolism and production that have equal values for the intrinsic rate of 
population increase, r. The arrows indicate possible trajectories when the organism is exposed to 
à toxicant, 


physiological state space. When the trajectory coincides with the local direction 
of the r isocline, there is no net effect on fitness (Figure 20.7). 

The theory of Sibly and Calow (1989) leads to the conclusion that the 
trade-off structure among life-history traits and physiological variables is of 
utmost importance for predicting effects of toxicants. Unfortunately, these 
trade-offs are not easily measured in practice, and so a test of the theory is 
difficult to make at the moment. There are many examples in the literature in 
which an expected negative correlation between life-history traits could not be 
demonstrated empirically (e.g, Van Dijk 1979; Hógstedt 1981; Meyer 1986). 
Ernsting and Isaaks (1991) showed that the costs of reproduction in the carabid 
beetle Notiophilus biguttatus depended on the way in which reproduction was 
manipulated (by temperature variation or by food availability). The analysis is 
further complicated by the fact that there may be differences between trade-offs 
at the physiological level and the genetic level. In the context of life-history 
evolution, it is the genetic correlation between two characters that matters. 
However, the existence of a negative genetic correlation between two traits does 
not preclude the existence of a positive phenotypic correlation between the same 
traits (Van Noordwijk and De Jong 1986). For the direct population effects of 
toxicants, the phenotypic trade-offs seem to be most relevant. However, for the 
selective pressure exerted by a toxicant, genetic correlations between life-history 
variables and toxicant tolerance are decisive (Posthuma et al. 1993). 


20.7 CONCLUSIONS 


The present chapter has emphasized the demographic approach towards 
ecotoxicity assessment at the population level. Other population approaches in 
ecotoxicology have not been addressed. Among these are the use of physiol- 
ogically-based population models (Hallam et al. 1990; Kooijman 1993; Klok 
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and de Roos 1996), the application of approaches developed for exploited popul- 
ations (Jensen 1984) and for populations of protected wildlife species (Grant et 
al. 1984; Emlen 1989), population analysis in multispecies test systems (Wolter- 
ing 1985; Ratte et al. 1992), and the incorporation of a spatial component into 
population level effects (Jepson and Thacker 1990; Sherratt and Jepson 1993). 

The material reviewed above has, nevertheless, illustrated the considerable 
progress recently made by the demographic approach. In several examples the 
study of population responses has provided a sharper picture of the effects of 
chemicals, although the looking glass effect mentioned in the introduction does 
not seem to have general validity. An important observation is that demo- 
graphic analysis may help to identify the ecological relevance of effects observed 
in individuals. The fact that the most sensitive variable is not necessarily the 
most relevant one in terms of r indicates the presence of plasticity and compen- 
sation mechanisms. These may be considered as the population analogue of the 
"functional redundancy" concept often discussed in an ecosystem context (see 
Levine 1989). The application of demographic analysis in an ecotoxicological 
context is a new area of research, and it opens up an interesting perspective for 
improving the ecological relevance of ecotoxicology. 
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